function [H_m1 H_m2] = func_H_m_l(q, m, CONSTS)

    k0 = CONSTS.k0;
    a = CONSTS.a;
    Q = k0*a*q;

    hankel_m_1 = hank_m_from_arg(Q, m, 1);
    hankel_mp1_1 = hank_m_from_arg(Q, m+1, 1);
    
    hankel_m_2 = hank_m_from_arg(Q, m, 2);
    hankel_mp1_2 = hank_m_from_arg(Q, m+1, 2);
    
    H_m1 = (hankel_mp1_1./(Q.*hankel_m_1)) -  (m./((Q).^2));
    
    H_m2 = (hankel_mp1_2./(Q.*hankel_m_2)) -  (m./((Q).^2));
    
end
